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In this paper we first consider a singular timelike spherical hypersurface embedded in a 
D-dimensional spherically symmetric bulk space-time. We analyze the different possibilities 
regarding the orientation of the gradient of the standard r coordinate in relation to the shell. 
The matter content of the shell is described by a number of non-interacting spherically 
symmetric matter fields. We analyze the dynamics according to Einstein's equations for 
matter satisfying certain energy conditions. Then we show that, under some circumstances, 
the system at a given proper time becomes unstable against an infinitesimal separation of 
non-interacting constituents. In particular, we consider the case where the shell is filled with 
collisionless particles that move in trajectories of constant angular momentum, and show 
that, for shells whose angular momentum distribution is discrete, the instability that we 
described before takes place. We also find this instability when the components are fluids 
with certain equations of state. We give explicit examples and construct solutions that 
represent a shell that splits into two shells. Then we extend those results for 5-dimensional 
Schwarzschild-AdS bulk space-times, which is a typical scenario for brane-world models, and 
show that the same kind of stability analysis and splitting solution can be constructed. We 
find that a widely proposed family of brane-world models are extremely unstable in this 
sense. Finally, we discuss possible interpretations of these features and their relation to the 
initial value problem with concentrated sources. 

PACS numbers: 

I. INTRODUCTION 

Singular surfaces proved to be a very useful tool in constructing toy-models to gain insight into 
general theoretical problems in general relativity, like the cosmic censorship conjecture [5]. They 
are also useful in constructing models for a number of concrete astrophysical and cosmological 
situations ([B], [7]) where shock-like surfaces formed, or where, as a valid simplifying assumption, 
some part of the system can be thought of as a singular surface. In this way, we consider that 
an analysis of general mathematical properties of this kind of objects is important and may have 
possible relevant astrophysical and cosmological implications. In particular, as a far as we know, an 
analysis on the non-uniqueness of the evolution of initial data associated with singular hypersurfaces 
has not been performed. 

Inspired by string theory, certain cosmological models that consider the visible universe as a 
part of a singular 4-dimensional surface embedded in a higher-dimensional space-time appeared 
[2], [3]. Besides being a theoretical possibility that is worth exploring, brane-world models have 
been proved to have interesting properties, such a Friedmann-like evolution equations that exhibit 
an inflationary stage They represent phenomenological models that can account for the "in- 
visibility" of extra dimensions, and constitute a framework where the holographic principle can be 
naturally understood. 

Outline 



We begin with a description of a general shell in a spherically symmetric bulk space-time of an 
arbitrary number of dimensions in Section|TTj In Section III we analyze some aspects of the dynamics 
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of these shells, provided their matter-energy content satisfy the dominant energy condition. In 
Section we develop the stability analysis against separation of non-interacting constituents and 
find a stability condition for general constituents, and then construct some examples of splitting 
solutions in Section |Vj In Section [VTt we analyze a splitting scenario that preserves Z2 symmetry in 
which a three component shell (with two of its constituents being identical) splits into three shells. 
Then in Section VII we extend both kind of instabilities to non-Z2-symmetric and Z2-symmetric 
brane-world models respectively. Finally, in Section VIII we comment on possible interpretations 
of these features in the context of the general initial value problem with concentrated sources. 
We formulate the following conjecture: unstable hypersurfaces are not thin-shell limits of thick 
configurations. 



II. SPHERICALLY SYMMETRIC THIN SHELLS IN D = n + 2 DIMENSIONS 

We begin with a description of a singular shell embedded in a spherically symmetric vacuum 
space-time of an arbitrary number of dimensions. The idea of spherical symmetry is originally 
defined for n + 1-dimensional riemannian manifolds: there is a group of isometries isomorphic to 
the SO{n) group. For lorentzian manifolds, we say that there exists a spacelike foliation such that 
every riemannian slice is spherically symmetric and the lapse and shift are invariant against the 
action of the group. Charting the space of group orbits with orthogonal coordinates (xq, xi), which 
is always possible, we may write 



-/(xo, xi)(ixo + /i(xo, xi)dx{ + r{xo,xif'dVL^. 



(1) 



We have a singular timelike orientable hypersurface S embedded in the bulk space-time. Because 
of the symmetry, the surface should be described by an equation $](x*^,x^) = 0. In particular, in 
gaussian coordinates, in a neighborhood of S the metric reads. 



ds^ 



-fiT,rj)dT^ + dr]^ + r{T,rjfdnl 



(2) 



where t/ = characterizes the surface, and /(r, 0) = 1 (r is the shell proper time). 

Einstein equations imply junction conditions on the surface [1] that relate the jump of the 
surface extrinsic curvature with the effective mass-energy tensor on the shell. These are the so- 
called Israel-Darmois junction conditions. The extrinsic curvature in gaussian coordinates reads. 
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Idf Idr 

2 dr] ^~ ' r dr] ^~ ^ ' 



1 dr 
r drj 



|r)=0 



(3) 



where, and from now on, latin indexes represent coordinates (r, 0i,..,0„) on the shell. In these 
coordinates the intrinsic metric takes the form. 



dsi 



-dr^ + Rirfda 



(4) 



where R{t) = r{T, 0). 

By virtue of Birkhoff theorem, vacuum solutions of Einstein equations can always be written in 
the form. 



ds^ = -F{r)dt^ + F{r)-^dr'^ + r'^da 



(5) 



where F{r) = 1 - 2M/r""\ and M is a parameter equivalent to the Misner-Sharp mass inside the 
innermost group orbit. From now on, we assume i2(r)"'~^ > 2M, unless otherwise stated, since 
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we are not interested in collapsing situations. With these expressions, we can write the extrinsic 
curvature on a given side of the shell in terms of R{t) and the parameter M that characterizes the 
bulk space-time there 



K'j = sign (^^\r,=o ] diag 



F'{R) + 2R ^R'' + F{R) ^R-^ + FjR) 
R R 



2JR^ + F{R) 



(6) 



In this way, giving M and specifying whether r increases or decreases with rj, we get an expression 
for the extrinsic curvature in terms of the intrinsic function R{t). 

On the other hand, the matter content of the shell is typically described by a tensor Sj defined 
on S such that we can formally write the Z)-dimensional stress-energy tensor as, 

= d{^)ss (7) 

where the S tensor in shell coordinates takes the form 

S} = diag[-p{T),p{T),...,p{T)]. (8) 

So, as a result of the symmetry imposed, we can describe the matter content of the shell as if 
it were composed of a n-dimensional perfect fluid, whose flow lines follow the trajectories of the 
comoving observers. Explicitly, 

s'^ = ph'^ + {p + pyu^ (9) 



where hij is the intrinsic metric defined in (|4j), and = {d/dry. If there is not hysteresis, we 
should be able to write p and p as functions of R. In that case, conservation of the source would 
read. 

This equation together with an equation of state f{p,p) = 0, provided it exists, determine p{R) 
and p{R). 

Now we can see that Israel-Darmois junction conditions should relate R{t), it first two deriva- 
tives, and the parameters M with the matter functions p{R) and p{R). Concentrated matter on 
the shell implies the discontinuity of the extrinsic curvature. Looking at ([g]), we realize that the 
jump should be ascribed to a difference between the mass parameters for the empty regions at both 
sides of the shell, which we call M/ and M//, and, eventually, to different signs for dr/dr] at both 
sides. In this paper we consider situations where at least one of the two Schwarzschild regions, say 
region /, is interior, that is, it can be described by an inequality r < R{t) in terms of the standard 
coordinates for that region. The reason for this choice, which only precludes the case where both 
regions are exterior, will become clear when we analyze the junction conditions. From now on, we 
choose the r] coordinate to decrease when going into region I, which implies dr / dr]\^^Q- > 0. The 
junction conditions then read. 



n 
R 



±JR^ + Fh-JR^ + Fi] = -Kp 



F'rr + 2R F'r + 2R u - 1 / / ■ ^ ^ / ■ ^ ^ ^ 

2^/R^ + Fii 2^Jr^ + Fi R ' " 



(11) 



(12) 
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where Fi = 1 — 2Mj/r"~^ and it = sign{dr /dr]\^^Q+). The situation not considered by these 
expressions, the case where / and // are both exierior regions, is not generally regarded as physical 
because it would imply a negative energy density p. In the context of brane cosmology, however, 
a brane tension is assumed, which implies an effective constant energy density on the shell, that 
could in principle be negative, because it is a non-dynamical entity. But, if that were the case, we 
would not recover standard cosmology as a low energy (at least since recombination) limit. So even 
in those speculative scenarios, a negative energy density type of matter-energy is hard to justify 
phenomenologically. In this way, we will only consider matter that satisfies the dominant energy 
condition, which is a prescription that includes an eventual positive brane tension. In particular, 
in the case of a exterior region //, a positive effective energy density p would imply Mjj > Mj. 



From ( 11 ) we can obtain an equation of motion that results independent of the sign it and reads 



V{R) 



Ml + Mil n\Mii - MiY K^p'^R^ 



8n2 



(13) 



But, in obtainig this equation, we have squared some quantities, so it may contain spurious solu- 



tions, that is, there could be solutions of (13) which are not solutions of (11) with a definite sign 
for dr/dr]\^=Q+. Nevertheless, it can be shown that every solution of (13) is always a solution of 



one of the two versions of (11) (with a specific choice for it). Therefore, we can construct solutions 
of the junction conditions by solving (13) and specifying a posteriori dr / dr]\^^Q+ accordingly^. For 
a given solution of the equation of motion, R(t), the specification that we should made is the 
following 



sign{dr/di]\^=o+) = sign{n\Fi - Fn) - n^p^R') = sign{2n\Mii - Mi) - K>^i2"+^). (14) 
It can be shown that this specification is consistent, because the roots of 2n'^{Mii — Mi) — 



n-l 



< 



K p{R) R'^^ are always in prohibited regions (y{R) > 0) or inside an event horizon R 
2M^i ii-j . In this way, for a given initial configuration of the shell (for an arbitrary choice of param- 
eters {Rq, p{Ro), Mil, Ml)), we can tell whether region II should be interior, or exterior, in order 
to satisfy the junction conditions, and that specification turns out to be consistent throughout the 
evolution ^. Now it is clear that giving a function p{R) (something that we can find for example 
imposing an equation of state) and the parameters {Mi,Mii), we have a complete description of 
the possible motions of the shell by solving an equation of motion for a single particle moving in 
a given potential. 



On the other hand, equation (12) implies. 



R 



n 



2R 



1 + R^T 



2n 



n 



1 



a + 1 A/l + i?2 



2Mi 



l + i?2 



IMii 



(15) 



where a = p/p. This expression will be useful later on. 



III. EQUATIONS OF MOTION AND MATTER MODELS 

In this Section we analyze general properties of the dynamics of the shell geometry, provided 
the matter content satisfies an energy condition. In particular, we analyze the dynamics for some 
matter-energy models of interest: fluids with equation of state p = ojp and Vlasov matter. 



^ The general problem of gluing two spherically symmetric space-times by means of an hypersurface was considered 
in [12]. 

Inside an event horizon, the r coordinate becomes timelike, space-time is not static, and the signature of dr/dr) 
may change without contradiction with Birkhoff 's theorem 
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Looking at the effective potential (13), we see that the permitted regions for the function R{t) 
are those in which V{R) < 0. The structure of these regions depend on the behavior of the 
function p{R), so we have to specify the matter model and its parameters in order to describe 
possible motions. The condition V{R) < can be explicitly written, 



(16) 



where it is clear that for R small enough, this inequality is always satisfied, so there are always 



"collapsing" solutions for the equation of motion (13) provided R{tq) is small enough. 

oo depends on the matter model as follows. 

then equation 



The asymptotic behavior of the potential when R 
If we suppose that matter satisfies the dominant energy condition, that is p > \p 



(29) implies 



So, giving an initial value po = p{Ro), we have for R > Rq 

Po{^y <p{R)<Po. (18) 

Considering the function a{R) = p{R)/p{R), the dominant energy condition then implies —1 < 
o:{R) < 1- We assume n > 2 and that a{R) has a limit Ooo when R goes to infinity, which implies 



that p acquires the asymptotic form p ~ CR^"^^^'^°°\ According to (14), we can see that if 
Ooo > —{n — l)/2n, for R large enough, region // would be exterior, and if Qqo < —{n — l)/2n 
then it would be interior. In the particular case Ooo = — ('^ — l)/^?^, region // character would be 
determined by the sign of 2n'^{Mjj — Mj) — k^C"^. The asymptotic behavior of V{R) in the limit 
i? — )• oo can be classified as follows (in all these cases A stands for different constants that depend 
on the parameters). 



• — 1 < Ooo < —{n — 

Then V{R) « -^^^^ ^ _^i^2(i-n(i+a^))^ ^^as case R{t) is unbounded, provided R{to) 
is large enough. 

• "oo = -(n - 1)/"- 

Then V{R) ~ ^" ^sn/ ^ ^ '^^^ condition for having an unbounded R{t) solution is 

an inequality involving the model parameters. 

• — (n — l)/n < Ooo < 

Then V{R) — )• |. There is always an upper bound for R{t). 

• Ooo = 

Then V{R) p ^ 2k? 2^2u~^''^ > cte. Once again the condition for having an unbounded 

R{t) is an inequality involving the model parameters. 

• < aoo < 1 

Then V{R) — " 2%"^!^^'' ^ —AR'^"'°'°°. There are always unbounded solutions. 
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Now it is clear that not for every matter model there exist mibounded solutions, so the limit 
i? — )• oo not always makes sense. Anyway, this limit, provided it exists, is important because of 
the fact that for large R we may find Newtonian or low-energy limits for the dynamics of the shell. 
For example in the context of brane- world cosmologies we hope to recover the standard Friedmann 
equations for R large in some sense. On the other hand, the existence of unbounded solutions will 
be useful to construct "unstable-against-separation" solutions, as we will see in the next Sections. 

It will also be useful to stress the fact that if we have k non-interacting matter fields (so we can 
write p = Yli=i Pi ^^'^ P — Yli=i Pi)i each with its own conservation equation and its corresponding 
ajoo, then Ooo = min{aioa}- In the next subsections we analyze some characteristics of V{R) for 
specific matter models. 



1. Fluid with equation of state p = ujp 



Equation (29) implies for this family of matter models 



p{R) = pA^) ■ (19) 



In this case V{R) can be written 



^ ' 2 2i?"-i 2C2 8n2i22(n(i+c.)-i) ^^"i 

where C = kpqRq^^^'^^ The dominant energy condition then implies — 1 < w < 1, and we have 
a{R) = Ooo = The character of region // should be determined by the sign of the expression 

2n^{Mii - M/)ii"(i+2<^)-^ - K^C'^. (21) 

We can see that this expression is a monotonically increasing function of R for u > —{n — l)/2n, 
a constant for uj = —{n — l)/2n, and a monotonically decreasing function for lo < —{n — l)/2n. In 
this way, for large R, in the first case region II should be exterior, and interior in the latter case. 
We now comment on some specific and relevant cases. 

• UJ = —1 represents a "cosmological constant fluid" or brane tension. We have p = po, and 
there are always unbounded R{t) solutions. For these unbounded solutions, region // must 
be an interior Schwarzschild region. For certain parameters {n, Mj, Mjj , kC) there are no 
prohibited regions. 

• UJ = represents dust, may be the simplest matter model of all. We have p = CR~^, and the 
condition for having unbounded R{t) is now Mjj — Mj > C^/n?, otherwise we would have 
an upper bound. Region // must be an exterior Schwarzschild region for the unbounded 
solutions. V{R) results a monotonically increasing function, so there could be at most one 
turning point, and if the parameters allow unbounded solutions, then there are no prohibited 
regions. 

• UJ = 1/n represents a radiation shell. We have p = CR^^^^^\ and there are always un- 
bounded solutions. For these unbounded solutions, region // must be exterior. Also in this 
case, for certain parameters there are no prohibited regions. 
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2. Vlasov matter 



Another simple and relevant matter model used for many astrophysical situations is Vlasov 
matter |15| . It has the advantage that its range of applicability is phenomenologically clear: a 
gravitating system can be thought of as an Einstein- Vlasov system whenever it is composed of 
many constituents with low collision probability. It is defined as an ensemble of autogravitating 
collisionless particles that interact with each other only by means of the curvature that they gener- 
ate as a whole. In particular, they follow geodesic trajectories. In this context, we have to specify 
what do we mean with "geodesic" because of the discontinuity of the metric connection on the 
shell. Furthermore, the curves contained within the shell can not be in general geodesies of neither 
connection (see [16j and \17\ for a discussion in the context of brane cosmology). The alterna- 
tive we have is to impose that the particles follow geodesies of the induced metric, something that 
sounds reasonable, specially taking into account that in this symmetry such geodesies are simply 
the trajectories of conserved angular momentum within the shell. The system we are considering 
in this subsection is a D-dimensional generalization of the four-dimensional case analyzed in ^ . 

If all particles are identical, the S tensor and the particle density current can be written, 

S'^ = -^, [ f{x,p)V^u'u^^^-^ , = [ f{x,p)V^u'^^^^ (22) 

J PO J PO 

where f{x,p) is the distribution of the number of particles in the tangent bundle and /i is the 
particles proper mass. Considering the symmetry and the conservation of the distribution of the 
number of particles with angular momentum modulus L (which we call n{L)), it can be shown 
that the independent components of the S tensor take the form, 

'■w = ^ / «(i)v'ra<iL , m = ^ / -l^dL (23) 

where Sn is the surface of a n-sphere of unit radius. 
Now we define, 

f{R) = j n{L)^J m + L'^dL (24) 

so the equation of motion for this shell can be written, 

T.rp^_l Mj + Mn n\Mjj - MifR^ Clf{Rf 

^ ' 2 2i?"-i 2Clf{RY 8n2i?2n ^^''> 

where C„ = nfi/Sn- 

Supposing that the function n{L) has compact support, then we have aoo = 0, and the asymp- 
totic behavior of V{R) would be 

R^O ViR)^-^^ (26) 
R^oo V{R) — ^ (27) 



n 



where (L) is the mean angular momentum modulus. We see that the shell can have unbounded 



motion if and only if CnN < n{Mjj — Mj). Considering derivatives of (25), it can be shown in 



particular that we can have oscillating or static solutions only if n = 2, which is the case analyzed 



in [9j. It can also be shown that (14) implies that region // must be an exterior Schwarzschild 



region for both the unbounded and the oscillating solutions. 
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IV. STABILITY OF A TWO-COMPONENT HYPERSURFACE AGAINST 
SEPARATION OF THE CONSTITUENTS 

We now consider a shell composed of two non-interacting matter fields. We impose that both 
matter fields are spherically symmetric, so their flow lines go with the comoving observers within 
the shell and they can be characterized as n-dimensional perfect fluids. The idea is to consider 
an infinitesimal separation of the components at a given proper time tq and calculate the relative 
acceleration of the separating shells. The resulting space-time, once both components are infinites- 
imally separated, can be determined merely by imposing continuity of Rq = R{to) and its first 
derivative, which is equivalent to impose smoothness (each "branch" should be an embedded sub- 
manifold) in the geometry of each "splitting" component, and also equivalent to impose continuity 
of normal vectors at the separation point. 

The stress-energy tensor takes the form, 

S'^ = {pi + P2W^ + {pi +P2+P1+ P2)u'u^ (28) 
Conservation of the source holds separately for each component, 

dj^ + <Pl±Pil=^ (29) 
Ctrl K 

SO we have in principle functions Pi{R) and Pi{R) for each component which are characteristic of 
the matter model. We should solve the equation of motion for the composed stress-energy tensor 



(28) and then for each component separately, taking into account the geometry that the new 
configuration would have, as illustrated in Fig.l. // the resulting relative acceleration turns out to 
have the same sign than the relative displacement, then the shell is unstable, otherwise it should he 
considered as stable. 

In this Section we prove that certain families of shells, that are solutions of Einstein equations, 
are unstable against this kind of separation, while other families are stable. What will remain to 
be proven is whether there are solutions that are initially stable and evolve as a single shell up to 
a point where they become unstable and consequently could undergo splitting. In the next Section 
we prove by giving examples that such situations are possible, and splitting shells that are formally 
solutions of Einstein equations can be constructed. 

The keys to analyze the relative acceleration of the resulting single-component shells are equa- 



tion (15) and Birkhoff's theorem. The space-time region that would appear between the splitting 
shells is characterized by a mass parameter Mint- We label with the index 2 {p2,P2) the splitting 
shell that moves into region //. The equation of motion for the shell labeled with the index 1 is 



obtained by replacing {Mi,Mjj, p{R)) in (13) with {Mj, Mint, pi{R))', and for the shell labeled with 
the index 2 by replacing those quantities with {Mint, Mjj , p2{R)) ■ Continuity of R{tq) determines 
Mint as a function of (Rq, Mj, Mjj) and the matter model parameters. 



Looking at (11), it can be seen that positivity of pi and p2 implies that Mint should be greater 
than at least one of the given mass parameters. In particular, it implies that if region II is exterior, 
then Mj < Mint < Mjj, and V^r in the intermediate region should point to region // (we label 
this possibility for the orientations with the A letter). In the same way, if region // is interior, then 
Mint > Ml J > Mj, and both orientations for in the intermediate region are possible (we call B 
the case where V^r points to //, and C the case where it points to /). In the latter case, in order to 



characterize the intermediate space-time, we must resort to the consistency condition ( 14 ) applied 
to one of the resulting shells. It can be shown that in any case the consistency conditions (14) 
applied for each splitting shell are indeed consistent with each other (in a given splitting scenario. 
Mint and V^r in the intermediate region are always well-defined). 



FIG. 1: Schematic representation of a splitting into two parts. The arrows inside the boxes represent the 
possible orientations of V°r in each region. We distinguish the three possible sets of orientations (which we 
call A, B and C) compatible with the positivity of pi and p2- 



As we said, we can distinguish three possible splitting scenarios, depending on the character 
of region II and the region that would form between the separating components, as illustrated in 
Fig.l. We want to compute for each situation the relative acceleration of the shells, and to achieve 
this we propose using (15) adapted to each shell as follows, 



Rl(ro) 



R2(ro) 



n-1 

" 2Ro 

n-1 

2Ro 



■2 f Sr 
1 + R (to) - s^gn 



-Qi + 1 



■ 2 / \ / dr 

1 + R (tq) - sign sign ] 

\ or]2 / \ or]2 



1 + R^(to)^ 



-02+1 



2Mj 



1 + R^{to)- 



l + fl2(ro) - 



l + i?2(ro)- 



2M„ 
Rn-' 



(30) 



where Oj = Pi/pi- We call the proper time coordinate of the original shell r and that of the 
resulting shells ri and T2, and make them coincide at the splitting point t = ti = T2 = tq. 

These derivatives are with respect to different time coordinates, so in order to make a comparison 
we should rewrite them in terms of a single time. A natural choice in a splitting scenario would 
be the standard time coordinate of the Schwarzschild intermediate region tint- To compute those 
derivatives we consider the following. For a general Schwarzschild space-time that has a shell as a 
boundary, and whose related variables and coefficients we note by a tilde, we can write. 



dt J F + R^ 



(31) 



where r is the proper time coordinate of the boundary surface, R its radius, and F = 1 — 2M/R^~^. 
This expression implies that in a splitting scenario we have dri/dtint = dT2/dtint at tq. In this way, 
at a given proper time tq we can write, 

^RM_U^^rY^^^^^^F^U U2F-r-,\ I (32) 



di'^ F \di J 2F 
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which in the sphtting scenario imphes 

We now can see that, as a consequence of the continuity of R{t) and it first derivative, the 
stability condition turns out to be simply 

sign (^Imt) (^2(ro) - i?i(ro)) < 0. (34) 



Taking into account (30) and (11) for each splitting shell, this inequality takes the form 

/3 = K\pi + p2)^i?"+^(n(ai + 02 + 1) - 1) - 2n^{Mjj - Mj){a2 - ai) > (35) 
which is independent of the signs of dr/dr]\int and dr/dr]\jj. It is remarkable that in all the three 



cases A, B and C, the stability condition takes the form (35). As a first result, we note that this 
condition implies that if both fluids are identical (a.2 = ai = a) and a > —{n — l)/2n, then the 
shell is always stable. A somehow unexpected result is that if we have two identical fluids with 
a < — (n — l)/2n, then the system is unstable. That would mean that if the matter content of a 
shell satisfies the latter inequality and it makes sense to separate the stress-energy tensor in two 
equal-alpha parts, then such a shell would always be unstable. We also note that if both q, > 
(which is the case for matter made of particles), then there is instability only if 02 is large enough 
in comparison to oi. In the next subsections we particularize the stability analysis for some specific 
matter models. 

Fluids with equation of state pi = ujiPi 
We propose two fluids with equation of state pi = ujiPi, where uj2 > <^i. In this case (35) takes 

> form 



the form 

(Cii?-"(^+'^i) + C72i2-"(i+'^^))2i?"+i(n(a;i + ^^2 + 1) - 1) - 2n^{Mii - Mj){uj2 - wi) > 0. (36) 

Provided wi > — (n — l)/2n, then the shell always become unstable for R large enough. 

As we have seen in the previous Section, at large R the oji fluid dominates p{R), and the 
asymptotic behavior is determined by Qqo = wi, so there are solutions with R arbitrarily large 
provided wi > 0. We are proving in this way that there exist shells so that, while being solutions 
of Einstein equations, are unstable and therefore, as we will argue later on, unphysical. 

Vlasov matter 

In the case of Vlasov matter, what characterizes the matter distribution is the function n{L). 
We may define non-interacting components merely by separating the particles into two sets that 
satisfy spherical symmetry on their own. For each set we would have an angular momentum 
distribution function nj(L) so that ni(-L) -|- n2{L) = n{L). We may define for each set, 

a- = 1 [ ML)L^ ,3^s 

' J ni{L)VWTL^dL J nVR^ + 

where we can see that < aj < 1/n. We commented that instability appears whenever 02 — ai is 
large enough. That could be achieved if one of the sets, say 77-2 (L), comprises all the particles with 
L greater than a certain value, as illustrated by the following case. 
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In particular, we analyze perhaps the simplest case that could display this kind of instability, a 



two-component shell where nj(L) = Ni6{L — Li). The condition (35) reads, 

Clf{Rf[{n + + nR^iLj + Lj) + (n - l)i2^] - 2n^R^+^{Mn - Mi){lI - L?) > (38) 

where we can see that for R large enough (A^ = A''i + 

• n = 2, the shell is stable, 

• n = 3, the shell is stable if K^f?N^ - 367r^(M// - Mi){Ll - Lj) > and unstable otherwise, 

• n > 3, the shell is unstable. 

As we discussed in the previous section, there are shells that go to infinity provided n{Mjj — 
Mj) > CnN. So a shell satisfying this last inequality and n > 3, or n = 3 and k^^j^'N'^ < 
367r^(M// — Mi){L2 — ), satisfies Einstein equations while being unstable at R large enough. In 
the next Section we give explicit examples where the shells are initially stable, become unstable at 
a later time, and where splitting solutions can be constructed. 

V. SPLITTING SHELLS 

In this Section we construct splitting solutions where the splitting takes place at a critical point 



of stability condition (35) (when the equality holds). These shells become unstable at that time, 
and a soft splitting solution may be built, that is, a splitting solution where R is also continuous 
(which is what defines the critical condition). 

Fluids with equation of state pi = wipi 

As a first example we consider a shell made of two non- inter acting fluids, whose parameters are 
as follows. 

• Example 1: n = 2, k = Svr, uji = 0, a;2 = 0.001, Ci = 0.01, C2 = 0.02, M/ = 1, Mu = 2 

We illustrate in Fig. 2 the velocity of the shell and the stability condition as functions of R. 
There are no prohibited regions for the motion of this shell, so if it is initially expanding then 
it will continue to expand forever (provided the system evolves as a single shell during its entire 
evolution). We can see that at Rc ~ 35,27 the shell becomes unstable. If we consider an initial 
radius i?o < Rc^ and say that the shell is initially expanding, then at R = Rc we may construct 
a splitting of the type illustrated in Fig.l, where the Mint parameter and the orientation of the r 



coordinate would be fixed by the continuity of R and the conditions ( 14 ) applied to each spitting 
shell respectively. As region // must be exterior, and the energy densities of each component must 
be positive, the r coordinate turns out to be a global coordinate (type C splitting). 

VJasov matter 

As a second example we consider a shell made of Vlasov matter, such that the particles only have 
two possible values for the modulus of their angular momentum {n{L) = Ni6{L—Li)+N2d{L—L2)). 
An example of this kind was already constructed in [9] for a 4-dimensional space-time (n = 2). In 
this paper we build a 5-dimensional example (like the previous one) whose parameters satisfy both 
n{Mij - Mi) > CnN and n'^^i^N^ < 367r^(M// - Mi){Ll - Lj). 
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FIG. 2: The solid curve represents as a function of R for a shell made of two non-interacting fluids with 
parameters n = 2, k = Stt, wi = 0, = 0.001, Ci = 0.01, Ca = 0.02, M/ = 1 and Mn = 2. The dotted 
curve represents 1000/3. The critical point is located at i? « 35, 27; for i? greater than that the shell results 
unstable. 



• Example 2: n = 3, k = 1, M/ = 1, Mu = 3, /xA^i = 30, /iiV2 = 23, Li = 2, L2 = 2.1 

We illustrate the effective potential V2{R) for this example in Fig. 3. There are no prohibited 
regions for this shell. We propose a separation betvi^een the particles with Li and the particles with 



-L2. At Rc ~ 27,98 there is a critical point of the stability condition (35), and for R greater than 
that the shell is unstable. As in the previous case, we propose a soft splitting at the critical point, 
provided Rq < Rc and that the shell is initially expanding, and obtain Mint from the continuity of 
R. In this case it also turns out that the r coordinate is a global coordinate. In Fig. 3 we illustrate 
the differences between the effective potentials associated with each splitting component and the 
one associated with the original shell. As Newtonian intuition suggest, we have that the shell 
composed of particles with angular momentum L2 begins to expand faster than the rate at which 
the original shell would expand, while the shell with Li begins to expand slower than that rate. 



VI. THREE-COMPONENT SHELL IN A Z2 SYMMETRIC SCENARIO 

In the context of brane cosmology, a common simplifying assumption coming from string theory 
is the existence of Z2 symmetry for the bulk space-time, such that the brane-world is the set of 



fixed points of this symmetry [3]. In this case, (14) confirms that region // must have the same 



character than region /, they both must be interior Schwarzschild regions with identical mass 
parameters. 

We consider a different kind of splitting in the Z2 symmetric scenario, as illustrated if Fig. 5. 
To preserve the symmetry we propose a fragmentation into three parts, one of which remains at 
the symmetry center, while the other two, being mirror images of one another, move away from 
that symmetry center. This is a three-component (non-interacting) shell splitting, where two of 
them are identical. We can consider the "outgoing" shells as two non-interacting identical parts 
of a single matter-energy component, whenever such interpretation makes sense. Analogously to 



Section IV, the key to analyze this kind of instability is to use (15) in order to obtain the relative 



acceleration of the resulting shells. We call M the mass parameter of the Schwarzschild regions 
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FIG. 3: V2{R) for n = 3, k = 1, M/ = 1, M// = 3, fiNi = 30, ^iiVz = 23, Li = 2, La = 2.1. There are no 
prohibited regions, and the curve is aheady very close to its asymptote at i? « 20. 
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FIG. 4: The dotted curve depicts V2{R; pi, Mj , Mint) — V2{R] p,Mi,Mii), while the dashed curve depicts 
V2{R]p,Mj,Mii) — V2{R; P2, Mint, Mil) . A solution representing a splitting like that illustrated in Fig.l 
can be constructed, where the separation takes place at i? w 27, 98. After the separation, the shell with 
L2 = 2.1 would expand faster than the shell with Li = 2. 



before the splitting, and M' the mass parameter that would correspond to the region between the 
symmetry center and one of the "outgoing" shells. We now define pQ and pQ as the stress-energy 
tensor components of the central splitting shell, and denote the corresponding components for the 
outgoing shells by p/2 and p/2, so we can write p = po + p and p = Po + p- The equation of 
motion of the resulting central shell is obtained by replacing {Mj, Mu , p{R)) with (M', M' , pq{R)) 
in (13), and for the outgoing shells by replacing those parameters with (M' , M, p{R)/2) in the 
same equation. Continuity of R{tq) determines M' as a function of (M, Rq) and the matter models 
parameters, and it can be shown that M' is always well-defined. Positivity of each energy density 
component implies that V^r always points to the central brane, which in turn implies M' > M. 
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FIG. 5: Schematic representation of a splitting into three parts that preserves Z2 symmetry. The arrows in 
each empty region represent the orientations of VqT. 
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Analogously to Section IV, the condition for having stability would be simply 

%o)-i?o(To) >0. 



After manipulating (39) and taking into account that p > and po > 0, we finally get 

(2na + n - l)(po + p) + (2nao + n - l)po > 0. 
This inequality can also be written as, 

2n(po + Po) + 2(nd - l)pQ + (2na + n - l)p > 



(39) 

(40) 
(41) 
(42) 



where it is clear that if a > 1/n, then the shell is always stable. We now suppose that on the 
original shell there is brane tension plus matter satisfying the dominant energy condition, so we 
may write po = b + pm, Po = —b + Pm, which illustrates the fact that the resulting central shell 
must also have brane tension^. In this case, we can rewrite the inequality (41) as follows 



2(n(am + a + 1) - l)pm + 2{na - 1)6 + {2na + n - l)p > 0. 



(43) 



We can see that if < a < 1/n, and < drn ^ ^/n, then the shell is unstable for b >> pm and 
b » p. Hence, if the shell is composed of non- interacting particles (Vlasov matter) plus some 



^ If we consider that the brane tension component is indivisible, it must remain at the symmetry center. 
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other matter field {pm,Pm) and a dominating brane tension, then it is typically unstable against 
separation of components. In that case, {p,p) would be determined by (p3|)), and we may separate 
the particle ensemble in two identical parts (ni(L) = n2{L) = n{L)/2) in order to construct a 
splitting of this kind. In typical brane-world models, b turns out to be much greater than any 
other energy density at least since nucleosynthesis, so it is relevant to ask whether this kind of 
instability also appears in brane cosmology scenarios, which is the subject of the next Section. 



VII. SPLITTING BRANE- WORLDS 

In this Section we extend the instability analysis that we made in previous Sections to actual 
brane-world scenarios. We consider SMS type brane- worlds |10j in which the five-dimensional bulk 
space-time has a spacial 4-slicing in which every slice possesses a kind of 3-dimensional isotropy: 
spherical, planar or hyperbolic symmetry. For such space-time, every timelike surface defined in the 
quotient manifold of the group orbits is a FRW submanifold. Bowcock, Charmousis and Gregory 
|llj have shown that this kind of manifold, provided they are solutions of vacuum Einstein equations 
with a cosmological constant, are always 5-dimensional Schwarzschild-Anti de Sitter space-times. 
This is the main property that allows us to extend the analysis that we have done for spherically 
symmetric space-times, as it plays the role that Birkhoff's theorem played in Sections |IV| and |VI[ 
In this way, every open vacuum region in the bulk must be Schwarzschild-AdS, which is a static 
space-time characterized by two parameters: a mass parameter and the cosmological constant A5. 
The metric in standard static coordinates for this kind of space-time reads, 

where Fir) = k — 2M/r^ -|- r^/£^, k is the curvature index (1,0, —1) and l"^ = — 6/A5. 

In this context, the brane-world is a FRW timelike boundary layer separating two different 
Schwarzschild-Anti de Sitter regions (/ and //). We consider for simplicity the same cosmological 
constant at both sides, but eventually different mass parameters. Like in the previous Sections, 
this singular surface can be described by a timelike curve in a quotient manifold, so a proper time 
coordinate r can be naturally defined, and the extrinsic curvature in (r, x, 4') coordinates it is also 
expressed by ^ . This coordinate system on the brane-world submanifold are precisely the standard 
coordinates for FRW space-times, where the scale factor is simply a(r) = r(r, r/ = 0) = R{t) {rj 
is the gaussian normal coordinate). As in the previous analysis, the matter-energy content of the 
shell is described by two functions p and p, as expressed by ([s]). Then, Israel-Darmois junction 
conditions applied to this situation could also be written as ( 11 ) and ( 12 ). Analogously, we suppose 
that at least one of the vacuum bulk regions is "interior" , which means that it can be described by 
an inequality r < R{t) in the static coordinate, and which is also a consequence of imposing p > 
matter fields on the brane^. 



ds^ = -F{r)dt^ + F{r)-^dr'^ + 



+ X^dn'' 



(44) 



In this way, as follows from (11), there is an equation of motion for the brane like that of a 



particle with unit mass and zero energy in a potential given by, 

V(m-^ Mi + Mjj R^ 9{Mn-Mjf K^p^R^ 

2 2R^ 2£2 2KV2i?6 72 
and that equation of motion can also be written as follows 

2_(Ry^^^p\ HMii-Mif Mi + Mu 1 k 

\R 36 K^p^R^ i?4 £2 ^2 • 1^0^ 



Provided R > 2max{Mi , Mn}. 
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In order to recover the classical Friedmann equations as a low energy limit, brane tension is needed 



+ b,p = —b + Pm)- In that case, (46) can be written [T3], 



C 



+ 



-Pm 1 + 



26 



k Mi + Mil A4 



(47) 



where G = K^b/ASir, A4 = {ti^¥ + 6A5)/12, and C= 3{Mii - MiY/ldTrGb. We notice that there 
is a "dark radiation" term and a term related to the "asymmetry" of the bulk space-time which 
may cause an inflation- like stage at early (pm > b) times^. Because we consider a uniform bulk 



cosmological constant, (14) would determine whether a given solution of the equation of motion is 



actually a solution of the junction conditions as in the spherically symmetric case. We can see that 
in the case of Z2 symmetry centered on the brane, the consistency condition always holds (both 
empty bulk regions are identical, hence interior); while for general non-Z2 symmetric brane- world 



scenarios, ( 14 ) would typically imply that both bulk vacuum regions must also be interior, because 
of the observational constraints on dark radiation (which are constraints on the (M/ -|- Mii)/R'^ 
term). 



A. Instability against separation of components in non-Z2-symmetric scenarios 



We now propose a brane- world made of two non-interacting homogeneous and isotropic matter- 
energy components, so we may write p = pi + P2 and p = Pi + P2- One of this components could 
always be the brane tension pi = —pi = b. In this subsection we consider a splitting like that of 
Fig.l, and impose continuity of the normal vectors at the separation point. We assume that the 
whole space-time possesses the symmetries we imposed above, so the intermediate region has to 
be Schwarzschild-AdS, with some (M, A5) parameters. In this paper we are not considering the 
possibility of having different cosmological constants on each open vacuum region, so the geometry 
of the intermediate region would be determined by the mass parameter Mint and the orientation 
of Vr on it. As in Section |IV[ both Mint and Vr are uniquely determined by the junction and 
consistency conditions applied to each splitting component. 



The junction condition on the brane (12) implies. 



R 



1 

R 



27?2 

A: + — ^ + ^2 ^ (3a + 1) 



R'^ + FiJR^ + Fii 



(48) 



which can be applied to each splitting component in order to determine the relative acceleration 



of the separating shells. Taking the same steps as in Section IV, the stability condition turns out 



to be precisely (34), which in this context can be written, 

54(M// - M/)(a2 - ai) - k"^ p'^ R"^ {3{ai + 02) + 2) < 



(49) 



regardless of the relative orientation of Vr in the intermediate region. Now, if we consider that the 
brane tension stress-energy should be conserved, one of the two splitting shells should have brane 
tension. In this way, we may write ai = {—b+pini)/{b + pim) which in the low-energy limit simply 
becomes ai ~ —1, so the stability condition can be written 



K^6^i?^(3a2 - 1) - 54(M// - M/)(a2 + 1) > 0. 



(50) 



Because there are very tight observational bounds on dark radiation, this stability condition for 
"realistic" brane- world scenarios turns out to be merely 02 > 1/3, which is a condition that 



^ Which nevertheless can not account actual inflation, as explained in [14] 
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matter made of particles does not satisfy. We may consider {p2,P2) as the matter-energy content 
of the universe, or some non-interacting part of it (hke the dark matter component to a good 
approximation), which is thought to have always satisfied a < 1/3, and indeed it practically 
satisfies q f« in the matter-dominated era. We then conclude that non-Z2-symmetric brane- 
world scenarios are typically unstable against separation of non-interacting components. 

Although the name of this subsection makes reference to non-Z2-symmetric scenarios, this kind 
of instability would also hold for an initially Z2 symmetric bulk space-time. Condition (50) makes 
perfect sense if Mj = M//, which is equivalent in this context to having Z2 symmetry. In that 
case, a two-component splitting would necessarily break the symmetry, so if we have a reason to 
impose it, as it may happen in brane- world models, such a splitting scenario should not be taken 
into account. 



B. Instability against separation of a three-component brane-world in a Z2-symmetric 



scenario 



Nevertheless, as we did in Section VI, we may propose a splitting of a shell made of three non- 
interacting components that preserves Z2 symmetry, as illustrated in Fig. 5. The mass parameter 
of the intermediate region M' would be determined by the junction conditions and the continuity 
of normal vectors at the separation point. By means of an analysis completely analogous to that 



of Section VI the stability condition turns out to be, 

(3ao + l)/9o + (3a + 1) + po) > (51) 

where the index denotes the component that remains to be at the symmetry center, while the tilde 
denotes quantities that correspond to the matter conforming the two symmetrically "outgoing" 
shells. The brane tension component must remain within the central brane, so we can write 
po = b + pom and po = —b + pom- In the low energy limit, the stability condition for this splitting 
scenario turns out to be a > 1/3, which is the same condition than that for a two-component 
splitting. As we argued, if it makes sense to separate some non-interacting component of the 
universe into two non- inter acting and identical parts, something that in particular always makes 
sense for dark matter, then splitting branes can be constructed, that are solutions of Einstein 
equations suggesting a particular kind of instability for brane-worlds, as we will discuss later on. 



C. A "realistic" example 

In this subsection we want to illustrate with an example how to construct a splitting solution 
in this context. We consider a SMS brane-world universe whose parameters are in agreement with 
recent observations, interpreted in the concordance model framework. We consider a brane-world 
universe whose matter content resembles that of the A — CDM universe, so we may write. 



Pm = Pc{nc + nb){l + zf + el{l + z)^ (52) 

where pc = '^Hq /8ttG is the critical density at present cosmological time, Qc and Qb ai"e the relative 
proportions of dark matter and baryonic matter respectively at present with respect to the critical 
density, and is the present radiation energy density. For simplicity, we consider the simplest 
choice of parameters {k = M = 0) for both empty bulk regions, which particularly implies non- 
existence of "dark radiation" and Z2 symmetry (each bulk region is an identical piece of AdS^). 
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The Friedmann equation for such a universe reads, 

(53) 



H"^ = Hi 



where = ^4/(3-ffo)' ^^"^ have ignored every term other than the radiation term in the 
quadratic in pm part of this expression. For smah z, this equation takes the standard form, 

H^ = Hl{{nc + ^b){i + zf + nA) (54) 

which makes exphcit the relation 0,c + + = 1- 

We propose a splitting in three parts like that of Section VI, where the {p,p) component is the 
dark matter component {p = pc^d^+z)^ ,p = 0) and the splitting takes place at present cosmological 
time (z = Q). It makes sense to separate two identical dark matter components, which would 
constitute the two "outgoing" shells, as cold dark matter is supposed to be made of non-interacting 
particles, hence it can be modeled as Vlasov matter and in this case (at cosmological scales) simply 
as dust. We should emphasize however that we do not expect this kind of construction to have 
a priori any degree of correspondence with the actual universe., They are solutions of Einstein 
equations that correspond to initial value data that can be obtained from a family of brane-world 
models analyzed in the literature. As actual solutions, we may consider the phenomenology that 
one would expect for observers living in a splitting brane-world, just to explore the theoretical 
possibility. In any case, with this construction we provide another example of the lack of uniqueness 
in the evolution of initial data that corresponds to concentrated sources. 

As we commented, brane-world models are typically unstable against separation of components. 



Replacing p with p/2 and setting (M/ = 0, M// = M') in (46), the equation of motion for any of 
the "out-going" dust shells can be written as follows. 



^2_87rG. ^ 2hedr\ vrG ,3 , A4 - 47rG6 

where we replaced k and A5 as functions of G and 6, and 87rGerfr./3 = M' / R'^ is the "dark radiation" 
term associated with the mass parameter M' of the intermediate region. Taking into account the 
fact that p «h sX the separation moment, we can write 



H'^ = Hi 



^dr{l+zf 1 + 



2pc 



(56) 



where Q^r = ^dr{Ro) / Pc and z = Rq/R — 1. In this way, all we need to determine the equation of 
motion of the resulting dust shells is the parameter ^dr, which can be found by continuity of the 
normal vectors at the separation point, that is, by solving it from H{z = 0) = H(z = 0). From 
(54) and (56) we get. 



02 



i(2J^b + J7c) + l 



so (56) can be written 

H^ = H 



^ + (l + z)2 + ^(l + z)^ + l^A-A 

2pc 2 J 2 2pc 



This expression can be simplified neglecting terms (6 >> pc), so we can write, 



^2 = Hi 



l + —z[l + 

Pc 



(57) 



(58) 



(59) 
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where it can be easily seen that the dust shells almost immediately rebound {H = at z —p^/b), 
which implies that the shells scale parameter begins to decrease until they collapse and form a pair 
of bulk black holes with mass parameter M' . 

On the other hand, the equation of motion of the central shell after the splitting can be computed 
from by writing pm = ^^b(l + z)^ and Mjj = Mj = M' , 



= Hi [Ofe(l + zf + Oc(l + + Oa] 

where the second term of the right hand side corresponds to the "dark radiation" term M' /R^ 
this way, the motion of the central shell is determined by (54) for z > and by (60) for z 



(60) 

In 
< 0. 

As we can see, the time derivative of the Hubble parameter would be discontinuous for the central 
shell, and the dynamical effect of the "disappearance" of the dark matter component would be as 
if the dark matter of the universe transformed into dark radiation, as the initial {z = 0) effective 
energy density of dark radiation is precisely the dark matter density at the time of the splitting. 



VIII. CONCLUDING REMARKS 

As far as we know, this kind of phenomena in the context of brane-world cosmologies was 
not previously considered in the literature. Nevertheless, there might be reasons in the context 
of string theory to a priori preclude any chance of having matter moving off the brane, as one 
may consider branes as fundamental objects by themselves, which appear in the definition of the 
action that describes the theory one is seeking to build. We basically affirm that these systems, 
in the framework of General Relativity, are unstable, or at least they entail non-uniqueness in the 
evolution, and the existence of those features may be an argument against the plausibility of the 
brane-world models considered in this work, and of certain families of singular hypersurface models 
in general. We remark that the critical point for any of the stability conditions ( 50|51 ) takes place 



at very early stages of the universe, long before nucleosynthesis. In this way, the brane results 
unstable almost during its entire evolution, and becomes increasingly unstable as time passes. A 
"soft" splitting, like those illustrated in Section |v| where R is continuous, can only take place 
at the critical point of the stability condition, so it could only happen in the inflationary era 
(taking into account experimental bounds on the brane tension based on deviations of Newtonian 
gravitation at submillimeter scales [18j), or even before. We may speculate in further investigations 
the phenomenological consequences of having this kind of splitting in the context of the very early 
universe. 

In any case, this problem should be related to the confinement of particles in brane-world models, 
as it happens in the case of Vlasov matter in 4-dimensional space-time |^. In Section 5 of that 
paper, the stability against single particle evaporation was considered. It was then concluded that 
for two-component Vlasov shells (like the one described in the second example of Section |V]) the 
stability condition for separation of components is always stronger than the stability condition for 
single particle evaporation. We may conjecture that for Vlasov matter the confinement of particles 
is a necessary condition for stability against separation of components. Nevertheless, the particle 
evaporation analysis may be relevant on its own, as the concept of "particle" in the cosmological 
context may have more sense physically than the concept of "non-interacting-component" . Particles 
in the universe, both if they are dark or baryonic matter, may follow geodesies of the induced metric 
to a good approximation (as long as collisions can be neglected). But, as Anderson and Tavakol 
analyzed in [16] and we pointed out in [9], geodesies of the induced metric do not generally coincide 
with geodesies of the bulk metric. Indeed, the geodesies of the bulk metric on the brane are not 
uniquely defined because of the discontinuity of the affine connection: for a given initial point on 
the singular surface and a given initial tangent velocity there might be two geodesies (each one 



20 



going into opposed bulk regions), a single geodesic, or no geodesic at all. Gravitational confinement 
is related to the non-existence of bulk geodesic for freely falling observers on the brane, as it means 
that according to the connection of region II the geodesic goes into region / and viceversa^ . // there 
exist a hulk geodesic and there are no non- gravitational forces acting on that particle, then it should 
classically follow that geodesic. As analyzed in [17], the confinement of freely- falling observers to 
domain walls in Z2-symmetric models can be accomplished if and only if the matter content of 
the shell satisfy the strong energy condition, which is a condition that models with brane tension 
generally do not satisfy, so it is no surprising that these models turn out to be unstable against 
separation of components. 

We could also interpret the instability result as an evidence for the existence of some mechanism 
that would "prevent" the splitting, that is some interaction among the constituents that may 
compensate the instability. We may take into account for example a scalar field on the bulk, 
which is a widely considered scenario as it may naturally exhibit an inflationary stage and it is 
easily incorporated by string theoretical considerations. Another possibility is to consider Lovelock 
gravity for the bulk. These are possible future research directions, and it may happen that some 
brane-world scenarios turn out to be stable in this sense, so that the stability against separation of 
constituents could work as an argument in favor of some of them. As a first step, as we suggested, 
we may look into brane-world models that naturally incorporate particle confinement. 

Another interesting issue regarding these solutions is the lack of uniqueness in the evolution 
of boundary layers. We have seen different situations involving singular hypersurfaces where a 
splitting solution can be constructed. But the splitting prescription is something that we can not see 
by looking at initial data in Cauchy surfaces that intersect the singular surface before the separation 
point. For every splitting solution that we constructed in this work, there is a corresponding non- 
splitting solution where the shell moves according to a equation of motion, and it may collapse, 
oscillate or expand always as a whole. That initial data is exactly the same for the splitting solution 
as for the corresponding non-splitting solution, that is, we do not have uniqueness in the evolution. 
Furthermore, if we consider particle evaporation we could also built solutions where an atmosphere 
forms, creating even more possible evolutions. This result may not be surprising taking into account 
the distributional data that a thin shell represents. In other contexts involving hyperbolic PDE, 
like ideal fluids in Minkowski space-time, surface of discontinuity may form and evolve (shock- 
waves) and their evolution according to the original PDE may not be unique [19]. There is an 
external principle (an entropy law) that recovers uniqueness. In this context we may have the 
same situation, as one may identify the "physical" evolution among the mathematical possibilities 
by means of some principle arising from the specific physics of the matter-energy content of the thin 
shell. After all, the information about the splitting may be exclusively contained in the degrees 
of freedom associated with the constituents of the singular surface. A specific treatment of these 
issues should be considered in a separate paper. 

On the other hand, looking at the results in [20], where we shown that certain singular surfaces 
are indeed thin shell limits of thick configurations, and that all those singular surfaces are stable 
against separation of constituents, we make the following conjecture: unstable shells are not thin 
shell limits of thick configurations. If we consider thin shells as idealizations of ultimately thick 
configurations, those unstable solutions could be interpreted as bad idealizations of thick configura- 
tions, and hence they can only be taken into account if we consider the system to be fundamentally 
thin, that is, if there is some fundamental principle that simply precludes the possibility of having 
matter-energy moving off the surface. An unstable thin shell may be considered as a spurious 
solution arising from the lost of information associated with the fact that we neglected a spa- 



^ The case where the bulk geodesies are contained within the brane are typically of zero-measure. 
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cial dimension. A general demonstration of this conjecture constitutes a possible future research 
direction, and it should be deeply related to the general initial value problem for concentrated 
sources. In any case, the stability-against-separation analysis that we exposed in this work can 
be regarded as a novel analysis to be taken into account in order to constraint the plausibility of 
models involving singular surfaces. 
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